Laughlin spin liquid states on lattices obtained from conformal field theory 
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We propose a set of spin system wavefunctions that are very similar to lattice versions of the 
Laughlin states. The wavefunctions are conformal blocks of conformal field theories, and for filling 
factor v = 1/2 we provide a parent Hamiltonian, which is valid for any even number of spins and is 
at the same time a 2D generalization of the Haldane-Shastry model. We also demonstrate that the 
Kalmeyer-Laughlin state is reproduced as a particular case of this model. Finally, we discuss various 
properties of the spin states and point out several analogies to known results for the Laughlin states. 
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Phenomena in strongly correlated systems are gener- 
ally hard to understand and describe, and therefore sim- 
ple model systems exhibiting various behaviors are im- 
portant guides. Laughlin's wavefunctions [l| 
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e.g., have played a key role in explaining the fractional 
quantum Hall effect, and this has triggered interest in 
understanding the nature of these states in depth. (The 
complex numbers Z n are the positions of the particles 
in the complex plane, v is the Landau level filling fac- 
tor, and the product of the single-particle phase factors 
x(Z n ) is the gauge factor.) Finding parent Hamiltoni- 
ans for wavefunctions is also very useful because it tells 
us how the behavior described by the wavefunction can 
be associated to the interactions between the particles 
in the system, and it can guide us to find experimen- 
tal situations where such a behavior occurs, even if the 
Hamiltonian itself cannot be implemented directly. 

In 1987, Kalmeyer and Laughlin (KL) introduced the 
v = 1/2 bosonic Laughlin state on a square lattice Q. 
Such models are expected to play a crucial role in un- 
derstanding topological phases of lattice systems, in very 
much the same way as Laughlin states do in bulk. Fur- 
thermore, they may open up the door [3] for the ex- 
perimental realization and investigation of Laughlin-likc 
states under very well-controlled conditions in optical lat- 
tices. The KL state has been further investigated in 0,0] , 
but it has been a long-standing problem to find a par- 
ent Hamiltonian for the state. Within the last few years, 
Hamiltonians have been found that are exact in the ther- 
modynamic limit [6|-|8|. In the present paper, we take a 
different approach, in which we propose to slightly mod- 
ify the Laughlin states. The modification enables us to 
find a relatively simple Hamiltonian, containing only two- 
and three-body interactions, which is exact also for finite 
systems and arbitrary lattices. We demonstrate that the 
modified states are very close to the original Laughlin 
lattice states for finite systems, while they are exactly 
the same in the limit of an infinite square lattice. 



As in the spirit of [9j , the wavefunctions we propose are 
chiral correlators of conformal blocks. The key element in 
the derivation of the Hamiltonian is to exploit that this 
structure allows us to apply rules from conformal field 
theory (CFT). Using the properties of null fields in CFT, 
we have recently derived [1 01 ] nonuniform and higher spin 
generalizations of the ID Haldane-Shastry (HS) model 
[ll|,[l2j], and in the present paper, we extend these results 
further to obtain a generalization of the HS model to 2D. 
Previous work on finding parent Hamiltonians for the KL 
state on an infinite lattice has in part been inspired by 
the original HS model. Here, we complete this idea by 
demonstrating that the HS state and the KL state are in 
fact two limiting cases of the same model. 

In addition, we characterize the most important physi- 
cal properties of the modified states, including correlation 
functions, topological entanglement entropies, and entan- 
glement spectra. It has been noted numerically that for 
some fractional quantum Hall states, the entanglement 
spectrum corresponds to the spectrum of the CFT that 
defines the wavefunction on the boundary 14] . Utilizing 
the particular structure of the proposed wavefunctions, 
we are here able to show analytically that the entangle- 
ment spectrum for a two-legged ladder (as defined in [151 ] ) 
exactly corresponds to the ID CFT with central charge 
c = 1 for v = 1/2, and we trace this back to the fact that 
the Yangian symmetry is inherited at the boundary. 

Wavefunction The wavefunctions we propose de- 
scribe the state of N spin 1/2 particles at positions 
Zi, . . . ,Zn in the complex plane, where N is even. They 
are chiral correlators of products of vertex operators 
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mal ordering, a is a positive parameter, s n = ±1 is twice 
the z-component of the ro'th spin, and f(z n ) is the field 
of a free mass-less boson, i.e. [la ]. 
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Here, 5 S = 1 for J2 n s n = and S 3 — otherwise, and the 
phase factors Xp,s p can be chosen at will, since the chiral 
correlator is only defined up to a phase. For a = 1/2, 
we shall always choose Xp,s p — exp[«7r(p — l)(s p + l)/2], 
since this ensures that © is a singlet [17j |. 

For a — 1/4, we note that @ can be writ- 
ten (up to an overall phase) as a Slater de- 
terminant of the single-particle wavefunctions 
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■ ,N/2, (We use the convention that Iln^m 
is the product over m and n, whereas Ilm(^n) i s the 
product only over m.) The states ipk(z n ) are not 
orthonormal, but can be made so without changing the 
Slater determinant. We can therefore regard ([2]) as the 
state of N/2 noninteracting fermions. This simplification 
enables us to use exact numerical computations rather 
than Monte Carlo simulations for a = 1/4 when we 
compute properties of the wavefunctions below. The 
v = 1 Laughlin state can also be written as a Slater 
determinant, and indeed we shall see in a moment that 
a = 1/4 corresponds to v = 1. 

Connection to the Laughlin states We next investi- 
gate the statement that @ is similar to lattice versions 
of the Laughlin states, which are in turn closely related 
to the continuous Laughlin states. We expect the corre- 
spondence to be approximately valid for all lattice config- 
urations for which the distribution of the lattice points is 
not too far from uniform and also if the complex plane is 
mapped into other geometries. We shall here consider the 
case of a (finite) square lattice in the complex plane since 
it is mathematically convenient and the case of an ap- 
proximately uniform distribution on the sphere because 
this geometry eliminates all boundaries. 

It has been speculated in 18[ that the KL state is 
proportional to the conformal block in ([2]) with a = 1/2 
on an infinite square lattice. Here, we proof explicitly in 
the supplement [13] for the case of a 2AI x 2M square 
lattice centered at the origin and with lattice constant 
b = \/ Aita that the ratio between ^ with appropriately 
chosen Xp,s p and ((TJ with v — (4a) -1 is given by 
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up to an irrelevant overall factor, where fM(z n ) = 

( z n/b)Y\^(^n)( l ~ Zn/zm)" 1 - In particular, the two 
wavefunctions coincide for M — » oo. In brief, we 
prove this result by transforming the spins into hard-core 
bosons by writing s n — 2q n — 1, q n £ {0, 1}. This allows 
us to express ([2]) in terms of fM{z n ). We then take the 
limit M — > oo, compute foo(z n ) by algebraic methods, 
and compare the result to (p}. We note that the cor- 
rect density of particles is obtained by scaling the lattice 
constant rather than by changing the filling factor of the 
lattice. As a consequence, the lattice filling factor, which 




FIG. 1. (a) Comparison between — ln(|/5(z)|) and 
— ln(|/oo (z)|) = Tv\z\ 2 /(2b 2 ) + constant. The points almost 
fall on a straight line with unit slope (solid line), (b) Plot of 
ri m (^„) pnm\ as a function of n for 100 spins on a sphere. 



is always 1/2, only coincides with the Landau level fill- 
ing factor v for a = 1/2. Figure [Ha) demonstrates that 
\Jm(z) \ s=s |/oo(2)| already for a 10 x 10 lattice, and thus 
there is a close relationship between (JlJ and ([2]) even for 
small systems. 

The model can be mapped from the complex plane to 
the unit sphere with polar angle 9 and azimuthal an- 
gle <p by the stereographic projection z = v/u, where 
u = cos(6»/2)e^/ 2 and v = sin(6»/2)e~^/ 2 . A small com- 
putation shows that Py = {viUj — UiVj)^ 1 and also the 
wavefunction ^ are invariant under SU{2) transforma- 
tions of the pair (u, v). Note that dij = 2\pij\ is the 
shortest distance dij — | — rij | between spin i and spin j, 
where n; = (sin(#i) cos(4>i), sin(6 l i) sin(</>i), cos(#i)) is the 
position vector of spin i. Writing s n = 2q n — 1 as before, 
we find that ^ is proportional to the Laughlin wavefunc- 
tion on the sphere Iln<m Pnm qm/u 03 with v = (4a) _1 
except for an extra factor of Jln^m Pnm" ■ Figure [2(b) 
shows that | Ilm(^n) Pnm\ varies only little with n for 
= 100, and so the correspondence between the pro- 
posed wavefunctions and the Laughlin states is again ap- 
proximately valid (note that the phase of Ilm(^n) Pnm" 
can be absorbed in Xn,s„)- Here, and in the following, 
we choose the distribution of the spins on the sphere by 
minimizing X)j<j d~j 2 numerically. 

Hamiltonian For a = 1/2, the vertex operators can 
be regarded as representations of spin 1/2 fields in the 
SU(2) 1 WZW model. Using properties of null fields in 
this model and the Ward identity, we derive [ItJ a set of 
positive semi-definite and Hermitian operators 
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i = 1,...,N, which annihilate the state ©. In ((H), 
w ij = 9( z i)/(zi — Zj) + h(Zi), where g and h are arbitrary 
functions of Zi, and Si = (Sf,Sf, Sf) is the spin operator 
of the ith. spin. It follows that ^ is the ground state of 
H % and thus also of H = ^ H t /4+(N + l) £\ . S, ■ Sj/6. 
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FIG. 2. (Color online) (a) Low-lying part of the energy spec- 
trum for a 2 x 3 (+) and a 4 x 3 ( x ) square lattice centered at 
the origin. (The horizontal axis shows the degeneracy.) (b) 
Averaged spin-spin correlation function for 200 spins on the 
sphere as a function of the distance between the spins. The 
error bars are of order a few times 10 -5 , and so the results 
are only converged for d < 1.2. The symbols encode the sign 
of the correlations (plus for positive and circle for negative). 
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We note that H reduces to the HS Hamiltonian [i"l, 12| 
when z n = exp(2irin/N) and j«y = 2zi/(zi — Zj) — 1, and 
the construction is hence a generalization of the HS model 
to 2D and nonuniform distributions of the spins. In 17j . 



we show that if the ground state is degenerate, then the 
additional ground states will not satisfy the Knizhnik- 
Zamolodchikov (KZ) equation The KZ equation is 
derived from the Sugawara construction that builds the 
Virasoro generators in terms of the Kac-Moody currents. 
So a degeneracy would indicate theories obeying the Kac- 
Moody but not the Virasoro algebra. It would be sur- 
prising if such theories exist, and we therefore expect the 
ground state to be unique. Exact diagonalization of H 
for small systems (see Fig. Ufa) for examples) also sug- 
gests uniqueness. 

On the sphere, we can obtain a Hamiltonian, which is 
invariant under SU (2) transformations of (u, v) by choos- 

ing H = J2i |uir 2 (fl i (1) +fl i (2) ), where fff ] (#f } ) is ® 
with Wij = l/(zi — Zj) (wij — Zi/(zi — Zj)). Finally, we 
note that a Hamiltonian for the case a = 1/4 can be con- 
structed by summing single-particle Hamiltonians, each 
of which is the identity minus the sum of the projections 
onto the orthonormalized single-particle states. 

Properties To further demonstrate the closeness be- 
tween ^ and the Laughlin states, we compute various 
properties of © in the following. We note that all the 
numerical results presented below except those related to 
entanglement spectra are independent of X P , Sp - 

Spin-spin correlation function Since the systems we 
consider are too big for exact numerical computations, 
we use the Metropolis Monte Carlo algorithm to compute 
the spin-spin correlation function 

/ Q z QZ\ _ ^s 1 ,...,s N S ^s 1 ,..., SN (z 1 ,...,Z N )\ 2 
4 L S1 ,..., SJV \i>si,..., SN {Zl,---,ZN)r 

For a = 1/2, the state is SU(2) invariant and so (SfS*) = 
S a b{SfSj) t a,b — x,y,z. Figure fj^b) shows the average 



FIG. 3. (Color online) (a) Inverse correlation length as a 
function of a for 100 spins on a sphere, (b-d) Averaged spin- 
spin correlation function as a function of the distance between 
the spins for (b) a = 0.24, (c) a = 0.25, and (d) a = 0.26. The 
error bars in (b) and (d) are of order 10 -5 for all points, while 
the results in (c) are exact. The solid curve in (c) is minus 
the expression for the correlation function of the continuous 
v = 1 Laughlin state on the sphere found in 22]. 



of the spin-spin correlation function for 200 spins on the 
sphere and a = 1/2 as a function of the distance between 
the spins. For a given d, the average is taken over all spin 
pairs for which the distance dij between the spins falls 
within the interval [d — e, d + e[, where e = 0.005. As 
for the Laughlin state with v = 1/2, we observe anti- 
ferromagnetic oscillations and exponential decay of the 
correlations. 

The dependence of the spin-spin correlation function 
on a is investigated in Fig. [3] Except for a close to 
0.25, we find that the correlator decays approximately as 
exp(— d/£), where d is the distance between the spins and 
£ is the correlation length. Numerical estimates of £ _1 
are shown in Fig. E^a). We also find that antiferromag- 
netic oscillations occur above a = 0.25, but not below. 
This observation is consistent with the lack of oscillations 
in the correlation function for the v = 1 Laughlin state 
and the presence of oscillations for filling factors below 
unity. It is also consistent with the conjecture that the 
transition occurs precisely at v = 1 2(J 21 1. The tran- 
sition is illustrated in plots (b-d). Finally, an analytical 



expression for the correlation function of the continuous 
v = 1 Laughlin state on the sphere has been found in 
[22l |. and the figure shows good agreement on intermedi- 
ate length scales. 

Entanglement entropy The possibility of having 
quasiparticles with fractional statistics is a very impor- 
tant aspect of the Laughlin states, and it is therefore 
very relevant to check whether @ also has nontrivial 
topological properties. Here, we study the entanglement 
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entropy (EE), since this allows us to extract the total 
quantum dimension D 23] . More precisely, if we divide 
a system into two subsystems A and B and the system 
is gapped, the EE has the form aP — 7 + ... 24J, where 



(a) 10 



a is a constant, P is (proportional to) the length of the 
boundary between A and B, —7 is called the topological 
entanglement entropy and fulfils 7 = ln(D) [24, 25|. and 
the ellipsis stands for terms that vanish for P — > 00. For 
abelian systems, D 2 is the number of different quasiparti- 
cles one can get by fusing the fundamental quasiparticlcs 
in the system, and so D = y/l/v for the Laughlin states 
with l/v GN@. 

The linear increase of the EE with P is confirmed 
in Fig. Ufa). We use here the Renyi entropy S L ' — 

— ln(Tr(p^)), where pa is the reduced density operator of 
the spins in region A, because this quantity can be com- 
puted efficiently by use of Monte Carlo methods [26|, |27 1 . 
The boundary between A and B is assumed to be the cir- 
cle 8 — 6l, where 9l is defined such that the area on the 
sphere with 8 < Ol relative to the complete area of the 
sphere equals L/N. The length of the boundary is then 
proportional to P = (L/N) 1 / 2 (\-L/N) 1 / 2 , and the steps 
in the EE appear due to the discreteness of the positions 
of the spins. For the Laughlin states, there are more 
ways, in which the system can naturally be divided into 
two parts, and previous studies H, 29| have computed 
entanglement entropies for both orbital partitioning, in 
which region A involves a subset of angular momentum 
eigenstates, and particle partitioning, in which a particu- 
lar subset of the spins comprise part A of the system. We 
note that such complications do not appear here, since 
the spins are fixed at specific positions. 

The topological entanglement entropy cannot be read 
off reliably from Fig. |4}a) because small errors in the 
linear term due to irregularities in the boundary easily 
become dominant. Instead, we use the method proposed 
in 24} , which eliminates the linear term by considering 
a linear combination of entanglement entropies of suit- 
ably chosen regions. The results in [24j are for the von 

Neumann entropy, but it has been shown in [3(| that 

(2) 

7 = ln(_D) also holds for S L when region A is a topo- 
logically trivial region. For a = 1/4, we find 7 = 0, and 
for a = 1/2, we get -7 = -0.341 ± 0.057 from a Monte 
Carlo simulation involving 160 spins on a sphere. Both 
results are consistent with the expected values (0 and 

— ln(2)/2 « —0.347, respectively) and with the results 
obtained for lattice models in [31(. 

Entanglement spectrum It has turned out (typically 
from numerical studies) that the low-lying part of the 
entanglement spectrum, defined as the eigenvalues of 

— \u(pa), is often related to some theory on the bound- 
ary of A [32j]. Utilizing the particular structure of the 
wavefunctions ([2|), we can here derive such a connection 
analytically for the case of a = 1/2 and N spins dis- 
tributed uniformly on two rings in the complex plane, 
i.e., the spins are at the positions exp[27r(2m ± X)/N\, 
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FIG. 4. (Color online) (a) Renyi entropy S£> = -ln(Tr(pi)) 
for 200 spins on a sphere for a = 1/2 (blue circles) and a = 
1/4 (green crosses) when region A is chosen to be the L spins 
that are closest to the north pole. The solid and dashed lines 
are linear fits, (b) The coefficients in (|6} for N = 12. 



where n = T2, . . . , TV/ 2. Specifically, we prove in the 
supplement [17j that pa (the reduced density operator 
of the inner ring) is invariant under Yangian transfor- 
mations, which means that — ln(p^) is a linear combina- 
tion of the invariants in the HS model. More precisely, 
we can write — bx(pA) as a linear combination of the 
identity Hq, the two-body operator H2 = 2^n^m^™ ' 
S m z n z m /zl m , Znm = z n - z m , the three-body operator 

B3 i ^^ n -£ m -£p ' (Sm X Sp)z n Z m Zp/ {z nrn Z rn pZp n ^j , 

and operators with higher-body interactions, which we 
write as H r [33|]. Considering these operators as nor- 
malized vectors \Hi) with inner product (Hi\Hj) = 
TriHiH^/iTriH?) Tr(Hj)) 1 / 2 , we can write 



I - Mp A )) = c Q \H Q ) + c 2 \H 2 ) + c 3 \H 3 ) + Cr\H r ). (6) 

The coefficients are given for N = 12 in Fig. [4jb) . We 
note that the above results do not generalize to the case 
of more than two rings. 

Conclusion We have proposed a set of spin system 
wavefunctions that are in many respects analogous to the 
Laughlin states. The proposed mapping between spin 
states and Laughlin states builds on CFT and demon- 
strates the usefulness of CFT as a tool to gain insight into 
many-body systems. We believe that similar mappings 
can be found also for other quantum Hall states. In par- 
ticular, one can obtain spin analogies of the Moore- Read 
state by gen eralizing the higher level ID spin models pro- 
posed in [10| to 2D, which can be done straightforwardly. 

The analogy between Laughlin states and spin states 
can be carried even further since the method proposed 
in Q to incorporate quasiholes by introducing additional 
conformal operators can also be used for the state ((2]) . As 
for the Laughlin states, one can interpret the square of 
the norm of the spin wavefunction as a particular charge 
distribution, and we expect this distribution to screen the 
quasiholes. It follows immediately from the construction 
that the analytic continuation properties of the wave- 
functions are the same as for the Laughlin states with 
quasiholes. It would be interesting to investigate these 
ideas further. 
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Supplemental material 



Singlet property of the wavefunction for a — 1/2 

In this section, we derive the choice of phase factors Xp.s that ensures that the wavefunction in Eq. (2) of the 
main text is a singlet for a = 1/2. The easy way to see this is to note that the vertex operators used to construct 
the wavefunction for a — 1/2 coincide with representations of the spin 1/2 fields of the SU(2) 1 WZW model, and 
therefore the wavefunction is SU(2) invariant if the phase factors are chosen appropriately. Since the phase factors 
do not depend on z n , it is sufficient to determine them for z n = exp(2mn / N) . In this case, the wavefunction reduces 
to the ground state of the Haldane-Shastry model, which is a singlet, and the phase factors can therefore be read off 
by comparison. 

One can also proof the singlet property by direct computation as follows. We define the total spin operator 
S = Y^iLi Si with components (S X ,S V , S z ) as well as the step up (+) and down (— ) operators S ± — S x ± iS v . The 5 S 
in Eq. (2) of the main paper trivially ensures that S z \ip) = 0, where \ip) = J2 Sl s N 1 Psi,...,sn( z 1i ■ ■■■> z n)\si, ■ ■ ■ , sat), 
and therefore the question is how one should choose the phases in order to achieve that 5 ,± |^) = 0. Noting that 

SkW = 2 £ ' l f } s 1 ,...,s k - 1 ,-s k ,s k+1 ,... : s N {zi,.-.,Z N )\s 1 ,...,S N ), (1) 

8l,...,SN 

S IW)=--^ £ s k ii Sl ,.... Sk _ 1 ^ Sk . Sk+1 ,..., SN {z 1 ,... 7 z N )\s 1 ,...,s N ), (2) 



we have 



1 N 

(s U ...,S N \S \tp) = ^X^ 1 ~ S k)^s u ...,s k - I ,s k ,s k+U ..., SN (Z1,...,Z N ) (3) 
k=l 

N N N 

= 2 Et 1 - s^k^^i-i)"- 1 n (** z -y snsk n n (*» - ^r- w2 , (4) 

k=l Xk,s k ri(#fc) P=l n<m 

where S s ^ — 1 for X)n(^fc) s « ~ s k — ^ an d 5 S k — otherwise. Since the factor (1 — Sf.) forces s& to be minus one, 
we conclude that (s\, . . . , sn\S~\iP) is zero unless the spin configuration (si, . . . , Sjv| contains N/2 — 1 spin ups and 
N/2 + 1 spin downs. Implicitly restricting the analysis to such configurations, we can drop the delta function in the 
following. The last two products in Q can be discarded since they do not depend on k and cannot be zero. The 
question is then whether 

n n _(#fc_)(^fe- -Zn-) 

is zero, where p± (p = n or k) is the subset of indices p for which s p = ±1 and F(k-) = (— 1) — 1 Xk_,i/Xk-,— l- The 
dependence on z n+ is contained in the polynomial II n — z n + ), ano - smce z n + are arbitrary and the number of 
factors in Y\ n+ {zt_ — z n+ ) is N/2 — 1, we conclude that §5§ is zero in general if and only if 

^> (M n Z \zk -z ) (6) 

is identically zero for all q € {0, 1, . . . , N/2 — 1}. Let us define the map t, which transforms the numbers 1,2,..., N/2 + 
1 into the indices of the N/2 + 1 down spins, i.e. {t(l), i(2), . . . , t(N/2 + 1)} is the set fc_. Multiplying @ by 
Y\ n _ <m _{z n _ — z m _), the expression turns into 

N/2+1 N/2+1 

^ £ nmzi^-iy- 1 n o%o ( 7 ) 

p—l n<m 

(n,m^p) 
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Comparing this to the Vandcrmondc determinant 
V{ai,a 2 , . ■ . ,a M ) = 





al ■ 


■ a M 


a\ 


a\ ■ 


a M 


M-l 




■ a M ~ Y 


1 







(-1) M(M - 1)/2 II 



(8) 



Kn<m<M 



we get 



JV/2+1 



Dq = (_l)W2-D/4 £ F(t(p))^ (p) (-l)P-V(* t(1) ,^^ (9) 



p=l 



= (_i)N(N/2-l)/4 



7° 
Z 1 



7° 

*t(2) 
Z 1 



't(l) 



6 t(2) 



F(t(l))< (1) F(t(2))< (2) ••• + l)>: (f+l) 

+ 1) 

,1 



t(f +1) 



**(f+l) 



(10) 



For F(t(p)) 7^ and general Zt( p ), this determinant is zero for all q G {0, 1, 2, ... , N/2 — 1} if and only if F(t(p)) = 1 
for all p. The conclusion is thus that in order to achieve (si, . . . , sjv|S ,_ |V') = for all spin configurations (si, . . . , sjv|, 
we need to choose Xfe,-i = ( — l) fe_1 Xfc,i for fc = 1,2, ... ,N, which means Xk,s k — exp(i : |(fc — l)(sfc + 1)) up to an 
overall irrelevant constant factor. For this choice of Xk,s k , the wavefunction transforms as \ip) — > (— l)^ 2 ^) when 
all the spins are rotated by 7r around the x-axis. The same transformation transforms into S T , and therefore it 
follows immediately that also S + \ip) = 0. 



Proof that the proposed wavefunctions coincide with the Laughlin states 
in the limit of an infinite square lattice 

We now proof that the wavefunctions in Eq. (2) of the main paper coincide with the Laughlin states in the case of 
an infinite square lattice with lattice constant b — \fAira. The first step, which is valid for general z n , is to transform 
the spins into hard-core bosons by writing s n = 2q n — 1, where q n € {0,1}. Dropping irrelevant g„-independent 
factors, this gives 



JV JV JV JV 



where t) q = 1 for J2n=i ^ n = and <5 q = otherwise. Since 



71=1 



N N N 

<j q n ( z « - ^m)- 2aq - = ^ q n (- z ™r 2aqn n ^ - v*m)~ 2aaB 

= S^f[ (-^)- 2Q(N/2 - 9 ™ ) n (1 - Zn/Zm)- 2 "**, 



m—1 



we get 



JV 



JV 



JV 



JV 



n=l 



n=l 



(11) 



(12) 
(13) 



(14) 



where fmizn) = (z n /b) rim(^ra)(l — z n/ z m) 1 and the label M = v^V /2 is introduced, because we will specialize to 
the case of a 2M x 2M lattice centered at the origin in the following, i.e. z n = b{n\ + in 2 ), where n = (ni,n 2 ) and 
m,n 2 = -M + \, -M + §,..., M- \. 
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To compute foa(z) on this square lattice, we note that 



/00O + &) = (z/b+1) lim 

M->-oo 



(2/6 + 1) lim 
M-yoo 



n 



-M+i<mi,m 2 <M-i 



1 - 



n 



mi +27712 ^z/b-\-l 
-M+i<mi,m 2 <M-i 



n 



mi + irri2 

mi + im2 
mi — 1 + im-i 

m\ + im,2 



1 - 



z/b 



n 



= — lim 

M->oo - 1 - - 1 - mi — 1 + !W2 

-M+i<mi,m 2 <M-i roi+imj^/i 

-M-i<mi<M-§ 
-M+i<m 2 <M-i 



mi — 1 + im2 

z/b 



1 - 



mi + im.2 



n 



M - \ + im 



= lim 

m->oo J -- L M + ^ — %m 

-M+i<m<M-i 2 -Af+|<m<M-i 



n 



1 - 



2/6 



M - 5 + im 



lim 

M->oo 



n 



M - i + im- 2/6 
x M + \ - im + 2/6 



-M+|<m<M-| 

From the Taylor expansion of ln(l + to), where w is a complex number, we have 



1 + w = cxp I — 



p=l 



and therefore 



!oo 
E [(^m - z/b - l/2) p - (z/6 + 1/2 - zrnf] 

-M+J<TO<M-A P =1 



M 



lim exp 



E 



(2g+l)M 2 «+ 1 

9=0 -M+i<m<M-i 



2 (im-z/b- l/2) 2q+1 



Consider 



(im - z/b - l/2) 2q+1 = (im)' Zq+1 - (2q + l)(imf q (z/b + 1/2) + 0(m Zq - L ) 



\2q+l 



2q 



foo(z) 

(15) 



(16) 



(17) 



(18) 



The sum of m 2<?+1 is zero, since m 2q+l is an odd function of m, and so only the term proportional to m 2q contributes 
in the limit M — > 00. Combining p7|) and (IT5|) and using the relations 



M-i 

E 



, 2 " 



we find 



C(z) = lim exp 

M— ¥00 



E 



M29+1 

2g + l 



(-1)« 



9=0 



(M 2q ) and V ^1 = I 
^ 2q + 1 4 



(19) 



9=0 _M+i<m<M-§ 



^ m 2 «(2z/6+l) 



= exp 



--(2^/6+1) 



From this and a similar computation, we conclude 



r 7r 1 r 7r 1 

foo(z + b) = exp ^--(22/6+ 1)J foo(z) and f 00 (z + ib) = exp |--(-2iz/6+ 1)J /oo(2). 

Solving these equations (remembering that Re(z)/6 and Im(,z)/6 are half integers) leads to 

foo(z) ocexp{-i7r[Re(z)/6- l]Im(z)/6} exp[-7r|z| 2 /(26 2 )]. 



(20) 



(21) 



(22) 
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Inserting this into (| 14[) and comparing to the Laughlin states in Eq. (1) of the main text, we conclude that the 
wavefunction in Eq. (2) of the main text is proportional to the Laughlin state with filling factor v = (4a) -1 and gauge 
factor rin=i(Xn,2Q n -i(— 1)^ 2q (™~ 1 )«" exp{— 2aq n in[Re(z n ) /b - 1] Im(z„)/6}) for M — > oo. We note that the difference 
between this result for a = 1/2 and the gauge factor in [J] is due to the fact that the lattice points are at half-integer 
rather than integer positions. 



Derivation of the Hamiltonian 



The main idea in the derivation of the parent Hamiltonian for the spin wavefunction with a = 1/2, which we 
presented recently in a more general setting in [l(|, is to use null fields in the SU (2)i WZW model. A null field 
Xsi(zi) has the property that 

' ' ' <t>^ 1 {z l ~i)Xs i {z l )(j) Si+1 {z l+1 ) . . .(j) SN {z N )) =0. (23) 
Here, we specifically consider the null fields [l(| 



[Xsi(Zi)]a = ^EE \ SabS Si s' ( ~ i^2£abc(Sf) StS ^j [J b _ ) (?;), 



(24) 



where a, 6, c can take the values x, y, z, e a bc is the Levi-Civita symbol, (S£) SiS > are the matrix elements of the operator 
Sf, and J b _ 1 is one of the modes of the chiral currents in the SU(2) 1 WZW model. The next step is to use the Ward 
identity 

(<p Sl (z X ) ■ ■ ■ S( _! {z i ^x){J h -X < t>s' i ){Zi)(t>a i+1 (Zi+l) ■■■<t>s N (z N )) 

N (S b ) s . s ' 

= Yl Y^-^-^sAzx) ■ ■ - (t)s[{zi) ■ ■ - ^s'^Zj) . . .4) SN {z N )) (25) 
to rewrite (|23l) into an expression containing the chiral correlator, which defines the spin wavefunction. This gives 

(26) 



j&i) 3 1 \ b,c j 



Due to the singlet property of the wavefunction, we also have 



\ E ( S J ~ 1 E ZobcStS) \W) = \\ St + i J2 £ ^S\ ] |V) = 0, (27) 
where we have used the relations 

Si S^ — -S a b + ^ E £abc Si aR d '^^Sabc^ab'c' = Sbb'Scc' — Sbc'Scb' (28) 

c a 

in the last step. Making an arbitrary linear combination — g{zi) x (Eq. ([26]) ) + h(zi) x (Eq. ([27]) ) of the two results, 
where g(zi) and h(zi) are functions of Zj, we therefore conclude that the operator 

2 N 

Ci = - w iJ ( s i + iS i x s i) ( 29 ) 

with = g(zi)/ (z% — Zj) + h(zi) annihilates the wavefunction, i.e., Ci\ijj) — 0. A hermitian and positive semi-definite 
Hamiltonian can then be defined as 

Hi = Cj ■ C t -. (30) 

The result in Eq. (4) of the main text now follows by combining and (|3U)l and applying (|2"8|) . Note that the 
Hamiltonian in Eq. (4) of the main text is the generalization to 2D of the ID Hamiltonian in Eq. (69) of [lOj. 
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Uniqueness of the state that is both a ground state and a solution to the Knizhnik-Zamolodchikov equation 

In order for a state to be a ground state of the Hamiltonian H, it must be annihilated individually by all the 
operators Cf and J^i This means that the state must fulfil ([26]) and be a singlet. The wavefunction in Eq. (2) of 
the main text with a = 1/2 furthermore satisfies the Knizhnik-Zamolodchikov equation 



We now show that there is only one state, which fulfils both (|26|) and (I3T[) and is a singlet. To do so, we first multiply 
([3T| by Sf and use (gB) to obtain 



Using (|26|) . this simplifies to 



dzi 3 j-f \ 2{z t - Zj) 



i(*0 v j; 



Taking a = z and utilizing sf = 1, it follows that 



-^-ip sl ,..., SN (zi,...,z N ) = ^2 f SlS ° ^ s u ...,s N {z\,...,z N ), (34) 

OZi £\Zi Zj) 

but this set of equations has the unique solution 

^s u ..., SN (zi, ...,z N ) = c(si, . . . , sat) JJ(zj - Zj) S * S3 ' /2 > (35) 

where c is a function, which does not depend on z n . The requirement that \ip) is a singlet uniquely determines 
c(si, . . . , Sjv) up to an irrelevant overall constant factor. Therefore the state is unique and given by Eq. (2) in the 
main text with a = 1/2. 



Invariance of pA under Yangian transformations 
From Eq. (|29|) . we know that the operator Q; = + , where 



Q« = £ Srf, Q\ 2) = i *a S * x S ^ (36) 

j&Q J'(7«) 

annihilates the spin wavefunction when a = 1/2. We consider here the case Wij — Izij (zi — zj) — 1 = + z j)/( z i ~ z j) 
and take the N spins to be uniformly distributed on R rings, i.e., the complex coordinates z, are given by 

Zm=u r C, C = exp(2^/F), r = l,2,...,R, n = 1, 2, . . . , P, (37) 

where u r is the radius of the rth ring and P is the number of spins in one ring. Note that 
p p > p _ i p 

Esr~^ U r C, n + U r iC, n tt-^ U r ( n n + U r > ■r— ( U r CJ l + U r > „ . , , ocA 

W rn ,r'n' = > — Z= T^7 = > = = Frr ' ' r ^ r ' 3 8) 



and 



/-n — n' i i PI > n P1-. PI P1-. 

y «, rnrn , = y i — ,±1 = y _a_ + y _J_ = _ y _L_ + y _L_ = . (39) 

n(^n') n(#n') S n=l S n=l S ' n=l S n=l S 
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Defining F rr = 0, we can write 

p p p B P 

n— 1 n^n' r'f/rjn'-li-l r'=l n' — 1 



and hence 



RE P R P. R. P R P 



E E E^ 1 )"-^ = E E (^V' E ^ E -EE *w = s. (41) 

r— 1 r" — 1 n— 1 r— 1 r" — 1 r' — 1 n' — 1 r' — 1 n' — l 

The latter operator is the total spin operator. Since is a singlet by construction, we have = 0. It follows 
from this and Q rn \ip) = that \ip) is also annihilated by the operator 

r r p r r I p p 

E E E( F_1 )-'Qr'n = EE i F ~ l )rr' 2A r , + * E E E W r> n ,r»n> $r'n X S r „ n , ] , (42) 
r— 1 r' — 1 n— 1 r—lr' — l \ r"(/r') li'-l 

where 

A r = — ^ ' w rn rn i S rn x S rn ' (43) 

is the Yangian of the rth ring. For the special case R = 2, there is only one term in the sum over r" in (|42|) . namely 
r" = r. Since F , w, and the cross product are all antisymmetric under exchange of the indices, the last term on 
the right hand side of (142|) is zero in this case, and we get the simplified expression 



2 2 p 

l(2) 



E E E^ 1 )""^'" = 2(^- 1 )i 2 A 2 + 2(F- 1 ) 21 A 1 = 2(F- 1 ) 12 (A 2 - Ax). (44) 



r—1 r' — l n—1 



We hence conclude that (A 2 — Ai)|-0) = 0, and therefore exp(i(A 2 — Ai)t)\tp) — \ip), where t is some real parameter. 
From this the invariance of pa (the reduced density operator of ring 1) under Yangian transformations follows: 

p A = Tr B (|^}(^|) = Tr s \e^ A2 - A ^\ij}(^\e~ l(A ^ Al)t ] = e -* Alt Tr B [e lA ^\^){^\e- lA ^\ e lAlt = e~ lAlt p A e lAlt . 

(45) 

This result does not generalize to the case of more than two rings since there are more terms in the sum over r" in 
(H2J) in this case. 



